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SUMMARY 


A  two-dimensional  Laser  Doppler  Anemometer  (LDA) ,  using 
polarization  to  discriminate  between  the  components,  has  been 
constructed  and  used  to  measure  flowfields  in  a  ducted  coaxial 
jet  system  with  a  centerbody.  The  flowfields  studied  included 
central  jet  only,  annular  jet  only,  and  combined  central  and 
annular  jets.  In  all  cases  the  central  jet  was  CO2  and  the  annu¬ 
lar  jet  was  air. 

The  small  diameter  of  the  central  jet  nozzle  relative  to  the 
centerbody  and  duct  should  result  in  the  central  jet  behaving  as 
a  free  jet  for  quite  a  distance  downstream  (~  100  nozzle  diame¬ 
ters  in  this  system).  Studies  of  the  jet  velocity  profile  devel¬ 
opment  have  been  made  and  compared  with  other  studies  in  the 
literature.  The  velocity  measurements  have  been  combined  with 
concentration  profiles  to  yield  entrainment  by  the  jet.  The 
values  obtained  differ  from  previously  published  results.  The 
consistency  checks  available  in  this  system  have  demonstrated  the 
reliability  of  these  measurements. 

The  annular- jet-only  measurements  have  been  used  to  compute 
the  streamline  locations  for  the  flow.  The  recirculation  vortex 
established  in  the  wake  of  the  centerbody  has  been  mapped  and  its 
center  located  in  both  position  and  strength.  Consistency  checks 
of  the  mass  flow  across  the  duct  have  demonstrated  the  quality  of 
the  measurements. 

The  coaxial  flow  measurements  have  looked  at  how  the  central 
jet  perturbs  the  reverse  flow  profiles  of  the  annular  jet.  It 
was  found  that  as  the  velocity  perturbation  was  dying  out,  the 
effect  on  the  velocity  fluctuations  became  most  marked,  with  the 
rms  velocity  fluctuation  increasing  by  almost  an  order  of 
magnitude  over  the  measurement  with  no  central  jet.  Axial 
velocity  profiles  were  measured  for  several  central  jet  flows  and 


the  results  compared  with  previous  measurements  usinc  a 
one-dimensional  LDA.  With  the  normalization  used  the  results 
agree  very  well,  even  though  different  noz2les  were  used  in  the 
two  measurements. 
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SECTION  1 
I  UTRODL'CT  IOM 


This  three-part,  final  report  iocanents  the  research  program 
performed  for  the  Air  Force  Wright  Aeronautical  Laboratories, 

Aero  Propulsion  Laboratory,  by  the  University  of  Dayton.  The 
research  had  two  overall  objectives:  (a)  providing  profile  data 
that  can  be  used  to  evaluate  combustor  and  fuel  combustion  models 
and  (b)  evaluating  the  performance  of  combustor  models  and  dif¬ 
ferent  diagnostic  techniques  in  various  combustion  environments. 

This  volume,  Part  I  of  the  report,  documents  the  design  and 
development  of  a  two-dimensional  laser  Doppler  anemometer  and  the 
experimental  data  collected.  The  analysis  and  modeling  tasks 
involving  the  numerical  flowfield  predictions  and  their  com¬ 
parisons  with  the  experimental  data  are  described  in  Part  II. 

Part  III  describes  the  design,  development  and  performance  of  a 
two-channel  time-resolved  laser  Raman  spectroscopy  system. 

Study  of  the  fluid  dynamics-combustion  interaction  has 
become  increasingly  important  in  the  design  of  new,  high- 
performance,  efficient  jet  engines.  Designing  a  combustor  to 
meet  the  performance  needs  of  the  engine,  while  operating  on 
lower  grade  fuels  and  generating  fewer  environmentally  harmful 
products,  has  become  a  serious  challenge.  Programs  are  now 
developing  numerical  computer  models  which  can  be  used  to  design 
and/or  analyze  combustors.  To  validate  the  performance  of  these 
models,  experimental  measurements  are  required  in  combustors  with 
well-defined  boundaries  and  controllable  inlet  parameters. 

A  program  has  been  underway  for  several  years  to  produce 
laser-based  diagnostic  tools  for  nonintrusively  monitoring  the 
combustion  processes  in  an  axisymmetric  research  bluff-body 
combustor  (see  Fig.  1),  located  at  the  Aero  Propulsion 


1 


Laboratory,  Wr iqht-Patterson  Air  Force  Base.  The  scope  of  the 
program  includes  the  development  of  instruments  for  point 
Tieasuremen ts  of  the  velocity  of  the  flowfield,  magor  soecies 
concentrations,  and  temperature.  A  goal  of  the  instrumentation 
is  to  take  data  at  a  sufficiently  'nigh  data  rate  to  follow  the 
fluctuations  in  the  relevant  flow  parameters. 

Global  studies  of  the  environment  use  schlieren,  shadow¬ 
graph,  and  high-speed  photography,  and  are  made  under  the  same 
operating  conditions  as  the  point  measurements,  to  gain  a  system 
understanding  of  the  fluid  processes.  Concurrent  global  and 
point  measurements,  as  well  as  conditionally  sampled  point 
measurements,  will  oe  required  in  the  future  to  advance  our 
understand  mg  and  to  provide  the  details  needed  for  refining  the 
computer  models. 

This  report  details  the  design  and  operation  of  the  most 
recent  version  of  the  Laser  Doppler  Anemometer  (LDA).  The 
instrument  measures  two  components  of  velocity  and  was  used  to 
profile  the  velocity  flowfield  in  the  combustor  shown  in 
Figure  1.  The  flowfield  entrance  geometry  is  that  of  coaxial 
gets  with  a  central  fuel  jet  issuing  from  the  face  of  a  center- 
body  and  annular  air  entering  around  the  periphery  of  the  center- 
body.  A  recirculation  zone  is  established  in  the  wake  of  the 
centerbody  which  extends  about  one  centerbody  diameter 
downstream,  for  the  flow  conditions  being  reported  here  (see 
Fig.  2).  The  first  window  station  provides  optical  access  from 
50  mm  upstream  to  250  mm  downstream  of  the  centerbody  face, 
extending  past  the  end  of  the  recirculation  zone.  The  next 
optical  station  begins  at  400  mm  and  extends  to  "00  mm 
downstream.  Relevant  dimensions  of  the  combustor  are: 
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Studies  were  made  of  the  flowfield  for  the  system  it 
isothermal  flow  under  three  distinct  operating  conditions  speci¬ 
fied  in  Table  1.  In  the  first  test  case  the  flow  approximates 
a  free  jet,  since  Rp  <<  Rq.  The  velocity  profiles  are  coupled 
with  concentration  measurements,  made  under  the  same  operating 
conditions  with  gas-sampling  probes,  to  provide  information  on 
entrainment.  The  second  test  set  provides  the  baseline  flowfield 
for  investigating  the  annular  jet-central  jet  interaction  and  the 
influences  of  combustion  on  the  flpwfield.  The  last  test  set  in 
this  series  examined  the  interaction  of  the  two  jets  and  the 
resulting  flowfield. 


TABLE  1 

OPERATING  CONDITIONS  FOR  THE  LDA  MEASUREMENTS 


TEST  SET 

FUEL  FLOW 
C02  (kg/hr) 

ANNULAR  FLOW 
Air  (<g/s) 

REYNOLDS  NUMBER 
Re  (x  lO*3) 

1 

2 

0 

11 

6 

0 

32 

16 

0 

83 

2 

0 

2 

481 

3 


6 

2 

16 


2 

2* 

2* 


SECTION  2 

LASER  DOPPLER  ANEMOMETER 


During  the  last  decade,  Laser  Doppler  Anenonetry  has  been 
applied  extensively  to  the  study  of  fluid  notion.  LDA  affords  a 
means  of  measuring  fluid  velocity  non intrus ively .  The  tem¬ 
perature  of  the  fluid  is  rarely  a  serious  concern  while  optical 
access  to  the  measurement  volume  is  available.  As  a  result,  LDA 
may  be  used  in  combusting  .ows  where  other  sensors  would  be 
destroyed.  LDA  measurements  can  be  made  from  a  remote  position 
and,  if  needed,  by  using  a  single-ended  optical  con f in ur a t ion 
(sender-receiver  at  the  same  location).  Thus,  it  is  possible  to 
shield  the  optical  system  from  any  harsh  environment  near  the 
measurement  site. 

The  theory  explaining  the  LDA  technique  has  been  presented 
by  many  authors  (see  for  example  Ref.  1)  to  whom  the  reader  is 
referred  for  details.  The  most  commonly  used  optical  con¬ 
figuration  is  the  dual-beam  system,  in  which  two  beams  of 
equal  intensity  are  focused  and  crossed  over  at  some  location  in 
space  (see  cig.  3).  A  scatterer  passing  tnrough  that  location 
will  scatter  light  from  each  beam,  each  component  being  Doppler 
shifted  from  the  original  frequency  by  an  amount  depending  upon 
the  velocity  of  the  scatterer  along  the  direction  of  the  optical 
beam.  L'sing  a  square-law  detector  (e.g.,  photomultiplier  tube) 
to  monitor  the  scattered  light  results  in  a  homodyne  mixing  of 
the  two  signals  with  a  measurable  beat  note  being  produced 
(typically  DC-30  MHz).  The  frequency  of  this  signal  is  directly 
proportional  to  the  velocity  of  the  scatterer  along  the  direction 
perpendicular  to  the  bisector  of  the  incident  beam  directions  and 
in  the  plane  containing  them.  The  geometry  shown  in  Figure  3  has 
a  significant  advantage  over  other  geometries  used  for  LDA;  the 
Doppler  frequency  generated  is  independent  of  the  position  of  the 
detector.  As  a  result,  the  detector  can  collect  scattered 
radiation  over  an  extended  angle  and  the  signals  will  add 
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constructively.  The  LDA  system  used  in  the  studies  reported  here 
will  be  described  in  detail  including  the  design  criteria,  the 
electronics  for  signal  handling,  the  system  checkouts,  etc.  The 
electronic  hardware  interfacing  the  LDA  to  the  computer,  the 
special  electronics  processing  and  control  circuits  designed  and 
built  for  this  system,  and  the  computer  software  were  all 
reported  in  depth  in  the  Technical  Manual.* 

The  LDA  system  used  in  these  studies  is  two-dimensional, 
with  a  three-beam  optical  configuration,  and  using  polarization 
to  separate  the  two  orthogonal  velocity  components  measured.  The 
optical  system  is  shown  in  Figure  4  and  has  been  described  in 
detail  by  Mag  ill,  et  al . -  The  system  is  a  composite  of  two 
independent  channels  identical  to  the  channel  shown  in  Figure  3. 

The  laser  source  is  a  5W  Argon-ion  laser  operated  at  438  nm. 
The  maximum  laser  power  is  approximately  1W  in  this  single  line, 
more  than  adequate  for  these  studies.  Typically  150  mW  of  laser 
power  was  used  in  the  isothermal  studies. 

A  dry  refractory  powder  is  dispersed  in  the  flow  system  to 
act  as  a  scattering  source  for  the  Doppler  signal.  The  ability 
of  the  particle  to  track  the  fluid  motion  depends  on  several 
factors  including  particle  diameter,  shape,  and  density.  The 
submicron  particles  used  here  were  computed^  to  track  fluc¬ 
tuations  in  the  fluid  flow,  at  1  kHz  frequency,  with  better  than 
99  percent  accuracy.  For  the  mixing  studies,  two  seeders  were 
used  to  seed  each  fluid  independently.  The  seeders  had  to  be 
adjusted  carefully  to  match  the  seeding  rates,  so  that  the 
measurements  in  the  mixing  region  would  not  unduly  weight  the 
influence  of  one  of  the  fluids  over  the  other. 


*  Andrews,  R.  J. ,  A.  J.  Lightman,  and  ?.  D.  Magill,  "Operation  and 
Maintenance  Manual  for  a  Two-Dimensional  Laser  Doppler 
Anemometer,"  UDR-TR-82-152  (1982). 
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LEGEND 

1.  POLARIZATION  ROTATOR 

2.  SPATIAL  FILTER  AND  3EAM  EXPANDER 


to  Differentiate  the  Orthogonal  Velocity  Componen is . 
The  luminosity  monitor  is  also  shown  (O3). 
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The  scattered  signal  was  collected  at  10°  from  the  forward 
direction.  The  available  windows  (300  mm  x  76  mm)  on  the  flow 
system  limit  the  detection  region  to  near  the  LDA  optic  axis  when 
scans  are  made  in  the  flowfield,  so  that  shadowing  does  not 
degrade  the  signal  so  much  that  it  becomes  unusable.  Forward 
scatter  was  chosen  because  it  maintains  the  polarization 
signature  needed  to  distinguish  the  velocity  components  and  it 
provides  greater  signal  strength.  The  most  common  difficulty 
with  using  polarization  discrimination  is  crosstalk  between  the 
channels,  depolarized  scattering  from  one  channel  being  received 
by  the  other  channel.  In  tests  using  one  channel  of  the  LDA 
system  alone,  no  crosstalk  measurements  were  observed  in  the 
other  channel. 

Two  reasons  motivated  the  selection  of  polarization 
discrimination  for  the  LDA.  First,  the  optical  system  is  more 
straightforward  and  easier  to  harden  for  use  in  the  high  acoustic 
level  and  large  temperature-gradient  environment  where  the  LDA 
operated.  Second,  it  was  understood  at  the  beginning  of  the 
program  that  3-D  LDA  measurements  would  be  required  eventually. 
Using  polarization  for  2-D  LDA  made  the  other  strong  Argon  laser 
line  available  for  use  in  obtaining  the  third  velocity  component. 
Employing  color  separation  for  all  components  necessitates  using 
either  a  weaker  laser  line  for  the  third  component,  usually  the 
line  at  476.5  nm  or  496.5  nm,^  or  a  frequency-mixing  scheme 
having  degraded  velocity  resolution.4  Results  of  research  using 
this  technique  have  not  demonstrated  any  great  degree  of  success. 

The  three-beam  configuration  produces  two  orthogonal  LDA 
systems.  The  beams,  all  parallel  to  the  optic  axis,  pass  through 
the  crossing  lens  at  the  locations  shown  in  Figure  5(a)  with  the 
polarizations  shown.  The  three  beams  cross  on  the  optic  axis  at 
the  focus  of  the  lens.  Placing  a  target  at  the  crossing  loca¬ 
tion  produces  the  orthogonal  fringe  system  shown  in  Figure  5(b), 
both  sets  of  fringes  at  45°  to  the  optic  axis.  The  U  and  L  beams 
are  orthogonally  polarized  and  do  not  interfere.  A  photograph  of 
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the  fringe  system  is  shown  in  Figure  6.  In  normal  use,  the 
fringe  spacing  is  measure-'!  using  the  graticule  that  can  he  seen  i 
"inure  '■> ,  rather  than  relying  on  co-otte-i  oarameters .  "or  the 
conical  svstem  use-!  —  lens  focal  lencth  "11  on,  laser  wave  Leon  oh 
4RF  nm. ,  three  heats  on  a  iiateter  circle — the  measure  1 

fringe  soacing  was  10  -  0.1  in. 

The  M  bean,  the  optical  bean  connon  to  both  Lt>A  systems,  is 
frequency  shifted  by  an  acous to-optic  modulator  to  produce  a  fre¬ 
quency  offset  that  eliminates  directional  ambiguity.  The  Braga 
cell  operates  at  40  MHz,  permitting  measurements  of  Topoler 
shifts  of  30  MHz  before  interferences  occur  with  the  next  laser 
mode  at  higher  frequency,  or  with  the  laser  discharge  noise  at 
lower  frequency.  The  Dopoler  frequency  range  of  30  MHz 
corresponds  to  a  velocity  range  of  300  m/s  in  the  firact ion.  of 
measurement.  Th is  ■  adequate ly  covers  the  expected  velocities  with 
a  sufficient  safety  margin. 

2.1  OPTICAL  DESIGN  CRITERIA 

'""he  ootical  system  had  to  be  designed  around  the  physical 
Limits  of  the  combustor.  The  combustor  duct  has  a  2  34  mm  inside- 
dianeter  and  the  LDA  system  had  to  scan  across  it.  ""he  closest 
ootics  on  either  side  of  the  combustor,  the  sending  and  receiving 
optics,  had  to  be  separated  by  at  least  1000  mm  so  that  they 
would  both  be  clear  of  the  duct  flanges  at  either  extreme  of  the 
scan.  Thus  the  sending  and  receiving  optics  had  to  have  focal 
lengths  of  at  least  500  mm  for  the  optical  design  shown. 

The  LDA  spot-size  diameter  was  set  according  to  the  velocity 
resolution  desired  and  the  expected  velocity  range.  To  obtain 
0.5%  accuracy  in  measuring  the  Doppler  period,  it  was  necessary 
to  measure  the  period  for  32  cycles  (at  the  higher  velocities  the 
Lower  signa 1-to-noise  ratio  limited  the  gating  accuracy  to  about 
1 /6  fringe^.  The  highest  expected  velocity  was  30Q  m/s  at  the 
exit  of  the  fuel  tube.  ''or  a  10-um  spacing  cringe  system  this 
would,  provide  a  30-MHz  DopoLer  shift.  3y  selecting  the  correct 
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3ragg  cell  "node  of  operation,  the  resulting  measured  fregue.ocy 
would  be  70  MHz  (shifting  in  the  reverse  direction  results  m  a 
10-MHz  signal  and  a  commensurate  degradation  m  optical  spatial 
resolution).  For  32  Doppler  cycles  to  be  collected  CO  MHz 
signal),  the  scatter  must  reside  m  the  LDA  spot  for  at  least 
0.45  ys;  thus  the  minimum  spot  diameter  would  oe  135  un.  In 
fact,  it  must  be  considerably  larger.  For  the  fringe  system  as 
shown  at  t  45°,  the  minimum  diameter  is  16  0  ym.  If  the  fringe 
system  were  rotated  453  (to  align  with  the  Y  and  Z  reference 
axis),  the  minimum  diameter  would  increase  to  240  urn,  dictated  by 
the  zero  velocity  in  the  Y  direction.  To  increase  the  crossing 
region  from  which  acceptable  signals  could  oe  collected,  it  was 
decided  to  increase  the  spot  dimension  from  the  160  um  minimum 
size.  The  final  choice  was  an  engineering  selection,  determined 
by  the  availablity  of  an  off-the-shelf  lens  having  near  diffrac¬ 
tion-limited  performance.  The  LDA  spot  diameter  used  was  300  ym 
and  the  lens  had  a  focal  length  of  711  mm. 

Several  key  parameters  describe  the  LDA  spot  which  will  be 
used  in  this  discussion.  For  conciseness  they  are  listed  here 
( refer  to  Fig.  3) : 


Spot  diameter 

dm  =  (4  X/t)  f/De-2 

(1) 

Spot  length 

lm  =  dm/tan  9/2 

(2) 

Fringe  spacing 

dc  =  X/2sin  9/2 

(3) 

The  size  of  the  laser  beam  was  ad]usted  by  selecting  the 
collimating  lens  of  the  spatial  filter-beam  expander  (#2, 

Fig.  4).  It  was  chosen  so  as  to  pass  through  the  5-mm  aperture 
of  the  3ragg  cell  without  creating  any  interference  fringes  in 
the  LDA  spot.  A  diameter  of  ~1.5  mm  satisfied  this  requirement. 
The  beam  remained  close  to  this  diameter  to  the  focussing  lens. 
Using  Eq .  (1),  this  yielded  a  spot  diameter  of  294  ym  in 
agreement  with  the  measured  value.  The  fringe  spacing  from 
Eq .  (3)  is  9.32  ym,  marginally  less  than  was  measured.  This 

probably  results  from  the  crossing  point  being  slightly  displaced 
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from  the  exact  focal  length.  The  length  of  the  snot,  from 
Eq .  (2),  computes  to  11.3  mm.  The  parabolic  mirror,  M  in 

Figure  4,  is  centered  10°  off-axis.  If  the  scattered  light 
collected  is  limited  to  a  small  region  around  10°,  me  spot 
length  would  be  determined  by  the  projection  of  the 
200- ura-d iame ter  optical  fiber  through  the  crossing  length 
(adjusted  for  the  50-ura  blur  circle  of  the  mirror).  The 
collected  spot  length  would  be  1.44  mm.  In  fact  the  mirror's 
collection  angle  extends  from*  3.3°  to  16°  and  the  length  of  the 
LDA  spot,  which  had  to  be  measured  experimentally,  is  3  mm  (see 
Fig  .  7)  . 

2.2  ELECTRONIC  SYSTEM 

The  scattered  radiation  is  transferred  via  an  optical 
fiber0  from  the  test  cell  to  the  control  room  30  m  distant.  The 
optical  fiber  is  immune  to  electrical  noise  and  its  bandwidth 
exceeds  the  requirements  of  the  opjtical  signals.  The  fiber  over¬ 
comes  difficulties  with  coaxial  cables,  such  as  high-frequency 
attenuation  and  spurious  signal  pickup.  The  signal  is  brought  up 
to  two  matched  photomultipliers  and  from  there  to  two  modified 
signal  processor  units  through  preamplifiers  (see  Fig.  3).  CPg 
computes  the  Doppler  period  for  a  preset  number  of  cycles,  with 
or  without  comparison,  as  selected  (see  TSI  1990A  manual)  and 
also  counts  the  number  of  Doppler  cycles  in  the  burst  (N-cycle) 
to  a  maximum  of  (Standard  1990A's  are  limited  to  2®-l 

cycles  and  could  overflow  once  velocities  dropped  below  47  m/s  in 
this  system,  losing  the  residence  time  information. )  CPg  com¬ 
putes  the  Doppler  period  only.  Both  processors  signal  that  data 
are  ready  on  their  data  bus  at  the  end  of  the  complete  Doppler 
burst.  If  the  data  rate  is  too  high  at  the  location  of  interest, 
the  processors  can  be  limited  in  their  data  acceptance  rate 
(unlimited,  30,000,  10,000,  3,000,  1,000,  300,  100  per  second), 
using  a  specialized  control  circuit  designed  and  constructed  by 
the  University  of  Dayton.  Another  detector  monitors  the 
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luminescence  along  the  line  of  the  anemometer's  optical  axis 
{used  in  combustion  studies  only).  T.nis  signal  is  processed  ay  a 
band-pass  filter  and  digitized.  The  digital  data  are  packed 
together  with  the  M-cycle  data,  4  bits  luminosity,  plus  12  cits 
M-cycle  data  in  tne  16-bit  word  fed  to  t.ne  computer. 

A  10-MHz  clock  measures  the  interarrival  time  of  the  signals 
on  each  channel  (Figures  3  and  9).  As  configured,  tne  clock  has 
a  12-bit  resolution  and  a  27-oit  dynamic  range,  permitting 
measurements  to  409.5  us  at  full  resolution  and  to  13.42  s  at 
reduced  resolution.  Figure  9  shows  the  digital  data  flow  on 
Channel  1,  where  three  16-bit  data  words  are  transferred  to  a 
fast  buffer  memory  and  by  direct  memory  access  (DMA)  to  tne  com¬ 
puter  memory.  Channel  2  transfers  the  elapsed  time  and  the 
Doppler  period  only.  The  DMA  unit  on  each  channel  can  process 
400k  16-bit  words/s,  permitting  a  long-term  data  collection  rate 
of  130k  data  points/s,  limited  by  Channel  1.  The  fast  FIFO 
memory  permits  storage  of  a  limited  number  of  data-points  when 
they  arrive  faster  than  the  DMA  can  process  them  (it  can 
accommodate  the  fastest  rate  of  the  counter-processor). 

2.3  DATA  COLLECTION 

The  sheer  volume  of  data  produced  by  the  LDA  system  necessi¬ 
tated  the  development  of  extensive  interactive  software  packages 
to  collect  and  analyze  the  data.  This  system  will  be  reviewed  in 
tne  following  paragraphs. 

The  LDA  data  are  collected  by  a  MODCOMP  minicomputer.  The 
internal  architecture  of  this  computer  restricted  the  data  to  one 
page  of  memory,  limiting  the  maximum  number  of  data  realizations 
to  about  5,000  per  data  point.  Typically,  this  provides  suf¬ 
ficient  data  in  tne  v,  w  coordinate  frame  to  ensure  the  stability 
of  the  computed  moments  of  the  distribution  (tne  velocity  com¬ 
ponents  u,  v,  w  align  in  the  X,  Y,  Z  directions  respectively). 
Once  the  data  are  collected  in  memory,  they  are  converted  from 
TSI-  and  CDRI -packed *  integer  formats  into  real  (floating-point) 


19 


V-'.M  . 


1 


& 


J 

•v 

f 


9.  Ir.e  Channel 
memory  13  r 
: CM A)  cent: 
shown . 


e  computer  Inter: 

.e.xed  th roach  e  itrec: 
7h  3  inter-  3  r  _  \  •/  =  1 


-  O  "1  _  o  - 


20 


numbers.  The  data  are  filtered,  analyzed,  reformatted,  and  put 
onto  disk  memory.  A  summary  of  the  data  (mean  velocities, 
standard  deviations,  data  rates,  etc.)  is  presented  to  the 
operator  on  the  computer  console.  In  addition,  a  visual  display 
of  the  velocity  histograms  and  a  summary  of  the  run  conditions 
and  results  are  plotted  on  a  Tektronix  graphics  terminal  (see 
Fig.  10).  If  the  data  are  to  be  retained,  a  hard  copy  is  made  of 
this  display  to  be  kept  for  reference.  The  average  velocities 
are  transformed  to  the  v  and  w  directions  and  are  passed  to 
another  computer  program  which  retains  all  data  sent  to  it  and 
can  plot  velocity  profiles  as  the  data  are  collected  (see 
Fig.  11) . 

When  the  on-line  data  collection  is  terminated,  the  data 
storage  disk  is  transferred  to  another  MODCOMP  computer,  referred 
to  as  the  HOST  computer.  The  data  are  further  filtered,  and 
pointers  are  attached  to  the  velocities  to  align  the  measurements 
of  the  velocity,  in  each  direction,  that  occurred  at  the  same 
instant.  (The  two  channels  operate  independently  and  not  all 
velocity  realizations  on  one  channel  necessarily  yield 
corresponding  realizations  on  the  other  channel.)  The  data  are 
then  reformatted  for  more  efficient  use  of  magnetic  tape  and 
stored  on  an  archival-data-record  magnetic-tape  set. 

2.4  DATA  PROCESSING 

Analyzing  the  data  obtained  by  LDA  measurements  is  a  subject 
for  active  discussion  by  researchers  in  the  field.  In 


*  To  speed  data  transfer,  the  information  has  been  restructured 
into  special  formats  that  occupy  less  than  half  the  number  of 
bytes  that  data  would  occupy  as  individual,  floating-point 
numbers . 
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high-volume  gas  flow  systems,  the  LDA  data  arrive  randomly. 
Counter-processors  interpret  the  Doppler  signals  and  their 
output,  which  is  discrete  in  time,  and  is  valid  only  during  t^e 
time  of  measurement.  Using  these  measurements  to  infer  t'oe  pro¬ 
perties  of  the  continuous  flowfield,  may  result  in  biased,  values. 
To  date,  a  variety  of  methods  have  been  proposed,  and  applied  to 
correct  biases  in  the  data  reduction  schemes,  but  the  underlying 
assumptions  of  these  methods  are  not  usually  tested  in  the 
experiment.  McLaughlin  and  Tiederman^  discussed  a  velocity¬ 
biasing  effect  in  measurements  of  flows  that  are  uniformly 
seeded.  Hoesel  and  Rodi7  considered  the  situation  where  the 
seeding  is  sufficiently  dense  that  the  LDA  measurements  occur 
rapidly  enough  to  follow  the  flow  accurately. 

Hoesel  and  Rodi7  and  Ruchhave,  et  al.^  considered  correc¬ 
tions  needed  when  the  time  between  realizations  was  longer  than 
the  correlation  length  in  the  flow.  Edwards9  pointed  out  the 
need  to  measure  the  integral  scale  length  at  every  location  to 
ensure  that  the  averages  presented  are  indicative  of  the  true 
mean  flow  and  not  just  local.  In  addition,  the  data  reduction 
and  interpretation  are  further  complicated  when  the  turbulent 
flow  has  a  time-dependent  mean  value. 

The  data  processing  programs  written  for  the  LDA  system 
incorporate  the  aforementioned  correctors  with  others  of  our  own 
design,  computing  the  mean  value,  standard  deviation,  skewness, 
and  kurtosis.  The  processing  schemes  include:  straight  arith¬ 
metic  processing;  using  a  McLaugh 1 i n-T iede rman ^  bias  corrector 
algorithm;  trapezoidal  integration  of  the  velocity  time  history; 
weighting  the  data  by  residence  time  ( Hoese 1-Rod i , 7  Buchhave®); 
simulating  processor-limited  data  collection  (supposedly 
resulting  in  truly  random  selection  of  the  data);l°  and  other 
schemes  to  reduce  the  possibility  of  overemphasizing  higher 
velocity  measurements.  Data  arrival  statistics  are  checked  to 
see  if  the  flow  is  uniformly  seeded.  A  test  is  made  of  the 
McLaugh 1 i n-T iede nan  assumption  by  analyzing  the  velocitv- 
interarrival-time  product  statistics. 


24 


▼ 


T 


T 


1 


r^T" 


J 


1 

i 

V* 

V. 

s'. 


The  processed  data  for  each  measurement  are  stored  on  a 
summary  computer  disk  (see  the  Technical  Manual  for  a  complete 
listing)  for  ease  when  searching.  At  present  the  results  of  all 
these  data  reduction  schemes  are  retained;  our  studies  have  not 
conclusively  validated  any  one  method.  Data  for  the  results 
shown  in  this  report  have  been  obtained  using  the  arithmetic 
averaging  scheme. 

2.5  SEEDERS 

Controlled  seeding  of  fluid  flow  for  LDA  measurements  is 
difficult.  Many  researchers  have  reported  seeders  which  operate 
well  in  their  restricted  environment,  but  the  requirements  of 
this  flow  system  have  proven  too  demanding  for  these  seeders. 

The  large  size  of  the  tunnel  which  needs  to  be  seeded  uniformly 
across  the  duct,  the  large  volume  of  air  and  seed,  and  the 
sustained  operation  for  long  periods  combine  to  make  demands 
beyond  the  capabilities  of  most  seeders.  In  addition,  the  flow 
system  has  a  coaxial  configuration  in  which  the  two  separate 
fluids,  central  jet  and  annular  jet,  must  be  seeded  to  the  same 
concentration  for  unbiased  LDA  measurements. 

A  commercial  seeder  (TSI  model  3400)  has  been  modified  for 
use  with  propane  as  well  as  other  gases,  and  our  operating 
procedure  allows  steady  operation  for  long  periods.  This 
procedure  includes  regularly  overhauling  the  unit,  using  only 
carefully  prepared  seed  (baked  dry,  ground,  and  sifted),  and 
regular  replacement  of  the  fluidized  bed  media.  In  addition,  the 
delivery  line  has  been  designed  to  supply  the  seed  within  the 
fuel  line,  away  from  walls.  Even  with  this  care,  the  range  of 
operation  of  the  unit  is  limited. 

Seeding  the  annular  flow,  which  can  exceed  2  kg/s,  is  more 
difficult.  The  seeder  must  put  out  10^  times  as  much  seed  per 
time  interval  to  balance  the  seeding  densities  in  the  annular  and 
central  flows.  A  high-volume,  fluidized  bed  seeder  that  uses  a 
cyclone  separator  to  remove  larger  agglomerates  has  been  built, 
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as  shown  in  Figure  12(a).  When  operated  with  a  solenoid 
the  air  supply  to  provide  an  air  hammer  stirring  effect,  this 
seeder  produces  sufficient  seed  to  balance  the  two  fluids.  To 
ensure  a  uniform  repeatable  cyclone  action,  the  fluidized  bed  is 
supported  on  a  wire  mesh  above  the  air  inlet.  An  unoalanced 
passage  mask,  shown  in  Figure  12(b),  forces  the  cyclone  action 
which  removes  the  larger  particles.  Another  smaller  diameter 
cyclone,  at  the  exit  port,  is  used  to  sift  the  seed  further. 

Then  the  seed  particles  are  dispersed  into  the  upstream  flow  in 
the  combustor  duct  via  a  dual-slotted  delivery  tube,  each  slot 
covers  30  percent  of  the  duct  diameter.  The  small  seeder  emits 
3  x  1()5  seed  particles  per  second  in  continuous  operation,  and 
the  large  seeder  emits  4  x  10^  particles  per  second  in  pulsed 
operation . 

Three  types  of  seed  material  have  been  tested  in  the  LDA 
system:  TiOg  with  a  flow  agent,  AI2O3  (alpha)  with  a  flow  agent, 

and  AI2O3  (gamma).  The  nominal  diameters  of  the  seed  particles 
are  respectively  0.2  urn,  0.3  um,  and  0.05  urn.  It  was  determined 
that  the  TiOg  seed  was  best  for  combusting  studies,  the  gamma- 
alumina  worked  best  in  cold  flow,  and  the  alpha-alumina  seed  made 
the  largest  contribution  to  depolarized  scattering,  probably  as  a 
result  of  its  jagged  shape.  As  mentioned  earlier,  Durst^  has 
shown  that  1.3  urn  diameter  TiOg  can  follow  1  kHz  fluctuations  in 
air  with  9  percent  accuracy,  and  0 . 4-u m-d iameter  particles  track 
10-kHz  fluctuations,  also  with  99  percent  accuracy.  Thus,  even 
if  the  Ti02  has  agglomerated  somewhat,  it  should  not  exceed  a 
size  which  would  distort  the  turbulent  energy  spectrum,  beyond 
the  measurement  uncertainties,  in  tne  frequency  range  of  interest 
( 1  to  5  kHz  maximum).  All  tne  seed  materials  should  provide 
accurate  spectra  over  the  entire  fre.;.ency  band  of  interest. 

2.6  LDA  SYSTEM  CHECK— LABORATORY 

The  LDA  system  was  tested  in  tne  laooratory  prior  to  being 
moved  into  the  tunnel.  The  first  parameter  validated  was  the 
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spot  diameter  (the  visual  checks  using  a  microscope  objective  to 
enlarge  the  spot  are  inaccurate).  A  chopper  run  by  a  synchronous 
motor  was  mounted  on  a  traversing  table  and  inserted  in  the  opti¬ 
cal  system  near  the  focal  position.  The  chopped  beam  was  picked 
up  by  a  solid-state  detector  (having  an  appropriate  bandwidth) 
and  displayed  on  an  oscilloscope.  As  the  chopper  blade  passes 
through  the  beam,  the  signal  will  go  from  off  to  on  (or  vice 
versa)  with  a  time  characteristic  that  depends  upon  the  laser 
mode.  For  a  TS.MQO  mode  and  a  chopper  blade  that  linearly 
traverses  the  laser  beam,  the  following  relations  are  valid: 

1  -  Ld)2 

Laser  Intensity  Profile  I  { r )  =  I^-^Je  a  (4; 

a  c 

Collected  Intensity  It(t)  =  I  F(vt/a)  (5) 

c  o 

where  IQ  is  the  laser  beam  intensity, 

F ( x )  is  the  error  function, 

v  is  the  linear  velocity  of  the  chopper  as  it  passes 
through  the  optical  beam, 
a  is  tne  beam  dimens  ion , 

and  t  is  the  time,  referenced  to  t=o  when  the  chopper  blade  is 
at  the  center  of  the  beam. 

Using  tabulated  values  of  the  error  function  the  10%  to  90% 
transition  of  the  transmitted  intensity  will  occur  in  a  distance 
of  2.56  a.  The  chopper  was  scanned  along  the  optic  axis  and  the 
minimum  rise  time  of  the  transmitted  beam  (Figure  13a) 
corresponds  to 

a  =■  76  um  t  4.5% 

The  commonly  accepted  definition  of  the  LDA  spot  diameter  is  the 
diameter  of  the  e-^  contour.  This  would  correspond  to 

dm  =  4a  --  30  4  um  t  4.5% 

as  designed.  By  traversing  the  chopper  along  the  optic  axis,  the 
beam  separation  and  hence  the  overlap  could  be  computed  from  the 


ooserved  oscilloscope  traces  (Figure  13b).  All  the  computed 
parameters  Hatched  the  design. 

Next,  the  LDA  signal  was  examined  using  a  scattering  source. 
The  3ragg  ceil  had  to  be  turned  off  to  allow  examination  of  the 
profile  and  the  Doppler  burst  together.  Typical  signals  are 
shown  in  Figures  14(a)  and  (b).  The  signals  have  the  textbook 
shape  and  quality,  showing  50%  visibility.  In  actual  operation, 
the  signal- to-noise  ratio  obtainable  is  much  less  than  for  these 
signals.  The  Bragg  cell  is  needed  and  the  bandwidth  of  the 
processor  must  be  expanded  to  encompass  the  expected  velocity 
fluctuations.  The  much  higher  frequency  results  in  a  decrease  in 
the  signal  in  direct  proportion  to  the  frequency,  as  a  result  of 
the  fewer  scattered  photons  collected  per  Doppler  cycle.  The 
expected  signal-to-noise  ratio  should  vary  as 

SNR  a  (f0Af)-1/2 

where  fQ  is  the  Doppler  frequency 
and  Af  is  the  amplifier  bandwidth. 

Scaling  the  SNR  observed  in  Figure  14  to  typical  operating 
conditions  will  result  in  a  minimum  degradation  in  SNR  of  15  to 
20  dB  (usually  more,  since  the  signal  shown  was  obtained  at  the 
center  of  the  LDA  spot,  the  optimum  location) .  Studying  the 
normally  obtained  LDA  signal  with  a  spectrum  analyzer  showed  that 
under  standard  operating  conditions  the  processors  operated  on 
signals  with  SNR  of  10  to  20  db,  sufficient  for  reliable 
process ing . 

2.7  LDA  SYSTEM  CHECK— TUNNEL 

The  LDA  system  was  tested  in  the  combustion  tunnel  before 
being  used  for  experimental  study.  Due  tc  time  limitations  a 
simple,  straightforward  test  was  desired.  The  simplest  test  was 
to  measure  the  fuel  exit  velocity  at  various  flow  rates  and  to 
compare  the  measurements  with  (a)  that  expected  from  a  simple, 
flat-top  profile  and  with  (b)  that  predicted  by  a  computation 
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that  took  temperature  and  compressibility  into  account.  The  fuel 
flow  was  increased  in  an  attempt  to  reach  choked  flow,  limited  by 
the  fuel  nozzle.  The  results  of  the  LDA  measurements  and  the  two 
predictions  are  shown  in  Figure  15.  In  the  region  M  <  3.5,  all 
three  are  in  excellent  agreement.  When  0.5  <M  <  0.67  t.ne  LDA  and 
corrected  predictions  agree.  3eyond  M  =  0.67  all  three  are  dif¬ 
ferent,  with  the  LDA  measurements  being  the  only  values  going  to 
the  correct  asymptotic  limit. 

The  quality  of  the  LDA  RF-bursts  observed,  the  close  match 
of  the  measured  optical  parameters  with  the  design  specifica¬ 
tions,  and  the  correct  low-velocity  measurements  leading  to  the 
expected  asymptotic  limit  (M=l)  all  gave  confidence  in  the 
operation  of  the  LDA  system.  After  these  checkouts,  the  system 
was  used  for  flowfield  investigations  as  described  in  the 
following  sections. 
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SECTION  3 
FREE  JET 


The  free,  turbulent,  axisym.metr ic  jet  is  a  standard  flow  in 
which  new  instrumentation  and/or  hypothesis  concerning  fluid  flow 
and  mixing  are  tested.  It  has  been  studied  extensively^  and 
there  exists  a  considerable  body  of  published  results  which  can 
be  used  in  any  new  study.  With  this  in  mind  it  was  decided  that 
the  free  net,  or  as  close  an  approximation,  as  can  be  obtained  in 
the  ducted  system  available,  should  be  the  initial  flowfield 
investigated  with  the  two-dimensional  LDA. 

The  central  ]et  issues  into  the  duct  through  a  nozzle  at  the 
center  of  the  bluff  body.  The  nozzle  diameter,  centerbodv 
diameter,  and  duct  diameter  are  in  the  relative  proportions  of 
1:29.3:53.2.  Thus  the  jet  will  have  the  characteristics  of  a 
free  jet  for  quite  some  distance  downstream.  The  short-taper 
nozzle  was  constructed  according  to  an  ASME  standard^  (see 
Fig.  16)  designed  to  produce  a  flat  profile  with  a  very  thin 
boundary  layer.  This  type  of  nozzle  was  selected  to  ensure  com¬ 
monality  with  the  work  of  other  researchers. 

Classical  studies  of  .the  free  jet  have  concentrated  on  such 
parameters  as  the  centerline  velocity  decay,  profiles  of  the  mean 
velocity  and  the  velocity  fluctuations  taken  at  a  fixed  axial 
location,  and  the  downstream  development  of  these  profiles. 
Recently  researchers  have  investigated  the  mixing  processes  in 
the  jet,  dividing  the  jet  into  different  regions  of  behavior  (see 
Fig.  17)  and  trying  to  interpret  the  parametric  variations  which 
characterized  these  regions. 1** 

The  study  of  jet  entrainment  has  recently  received  new 
experimental  attention  using  more  modern  tools.  Past  studies 
have  included  the  porous  wall  technique  used  by  Ricou  and 
Spaldingl5  and  Hill. 16  More  detailed  study  of  the  development 
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region  is  possible  by  examining  the  downstream  velocity  profiles. 
Crow  and  Champagne^-7  and  Boguslawski  and  Popiel33  used  hot-wire 
probes  to  measure  the  velocity  profile  in  the  region  out  to  10 
nozzle  diameters  downstream.  Their  results  have  been 
contrasted^  with  the  earlier  measurements  by  Hill  I'’  in  the  same 
region  and  the  further  downstream  results^-3  (out  to  several 
hundred  diameters).  Substantial  differences  have  been  found. 

The  LDA  system  was  used  to  study  the  velocity  profiles  in  the 
region  between  0.4  and  105  fuel  tube  diameters.  This  region 
overlaps  that  of  the  previous  studies  and  it  was  thought  that 
this  might  clarify  the  observed  differences. 

3.1  EXPERIMENTS 

The  experiments  concentrated  on  the  central  net  with  a  COg 
flow  of  6  kg/hr,  the  same  flow  used  in  the  coaxial  studies. 
Additional  profiles  were  made  at  selected  locations  for  CO2  flows 
of  2  and  16  kg/hr,  to  look  at  the  scaling  as  a  function  of  the 
central  jet  Reynolds  number  (see  Tables  2  and  3). 

TABLE  2 

FLOW  RATES  USED  IN  CENTRAL  JET  STUDIES 


Flow  (CO2) 

Reynolds  Number 

kg/hr 

Re  (xl0~3) 

2 

11 

6 

31 

16 

83 

'  - 
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TA3LE  3 

PROFILE  MEASUREMENT  LOCATIONS 


The  entrance  profile  was  approximated  by  the  profile  at 
Z/Op  =  0.4,  the  closest  approach  of  the  LDA .  The  central  jet 
boundary  layer  is  thin  in  this  region,  near  the  entrance  plane, 
and  the  velocity  gradient  is  steep.  The  LDA  spot  dimensions  form 
an  acceptance  aperture  which  is  convoluted  with  the  true 
flowfield  to  provide  the  measured  flowfield.  To  minimize  the 
effects  of  this  aperture  the  LDA  scans  were  made  in  the  vertical 
plane  going  through  the  centerline,  in  which  direction  the  LDA 
spot  aperture  is  minimum. 

The  measured  profiles  were  combined  with  density  profiles 
measured  with  intrusive  probes. ^9  a  numerical  integration  scheme 
(assuming  cylindrical  symmetry)  was  used  to  obtain  the  mass  and 
momentum  flow  within  each  profile.  The  mass  of  CO2  was 
computed  (the  intrusive  probe  measures  the  fraction  of  CO2  in  the 
sample,  the  concentration  being  diluted  by  the  air  that  the  jet 
entrains).  Mass  and  momentum  conservation  required  that  the  CO2 
mass  flow  be  constant,  equal  to  the  issuing  central  jet  CO2  flow, 
and  that  the  jet  momentum  be  constant  at  every  axial  station 
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(assuming  that  the  contribution  of  the  natural  abundance  of  CO 2 
in  the  entrained  air  is  negligible  compared  with  the  jet  flow). 

3 . 2  RESULTS 

Figures  13(a)  to  (m)  show  plots  of  the  central  jet  profiles 
measured.  Each  figure  presents  the  results  for  the  w  and  v 
components.  For  these  scans  w  is  the  axial  component  and  v  is 
the  radial  component,  except  for  13  (  i )  and  1 3  ( J< )  where  v  is  the 
azimuthal  component.  Each  page  has  four  plots:  average 
velocity;  rms  value  of  the  fluctuating  component;  skewness;  and 
kurtosis  of  the  measurement  distribution.  (Skewness  relates  to 
the  symmetry  of  the  distribution;  it  is  zero  for  a  symmetric 
distribution.  Kurtosis  is  a  shape  parameter  describing  the 
distribution;  it  is  3  for  a  Gaussian  distribution  and  1.3  for  a 
flat-top  distribution.) 

The  profile  at  Z/Dp  =  0.4  is  the  inlet  profile.  The  flat 
top  of  the  velocity  profile  is  evident.  Other  measurements  taken 
in  this  region  confirmed  the  flat  profile.  The  shear  layer 
thickness  observed  is  largely  due  to  the  width  of  the  LDA 
measurement  volume  (300  urn).  The  shear  layer  width  measured  from 
the  velocity  profile  indicates  the  true  shear  layer's  maximum 
width  is  on  the  order  of  100  urn.  The  symmetry  of  the  flow  (in 
the  measurement  plane)  is  evident  in  all  four  plots  shown. 

The  skewness,  when  traced  from  r/Rp  =  0  outward  shows  a 
sharp  drop  from  0  to  -1,  then  a  rapid  increase  to  +1,  and  then  a 
drop  back  to  0.  The  variation  occurs  within  the  observed  shear 
layer.  This  shape  is  anticipated  for  the  situation  where  the 
shear  layer  is  thinner  than  the  LDA  cross  section.  As  the 
measurement  volume  moves  into  the  shear  layer,  a  large  part  of  it 
remains  in  the  potential  core.  Hence  the  skewness  will  be 
negative.  When  the  spot  straddles  the  shear  layer,  the  skewness 
will  tend  to  zero.  As  the  spot  moves  to  the  outside  but  still 
overlaps  the  shear  layer,  the  bulk  of  the  LDA  volume  will  be  in 
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Figure  13. 


(a)  -  (m)  Profiles  of  Velocity  Distribution 
Descriptors  3ased  Upon  the  First  Four  Moments. 
Measurements  are  for  the  central- jet-cr.ly  condition 
with  CO2  mass  clow  of  o  kg  hr.  The  location  and  f'.ov 
conditions  for  each  set  are  given  at  the  ton  of  the 
page.  figures  _  3  f  1 )  a  no  1 3  t  <  )  are  X  —  scans  sc  the  ’ '  — 
component  is  the  azimuthal  component.  For  all  the 
other  plots  the  7-component  is  the  radial  component. 
The  W-component  is  always  the  axial  component'. 
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the  zero  velocity  region.  The  contribution  of  the  shear  flow 
will  provide  a  positive  skewness.  Once  the  spot  is  completely 
outside  the  flow,  the  skewness  should  decrease  back  to  zero. 

At  the  downstream  measurement  stations  the  shear  layer  crows 
in  width  and  the  effects  of  the  LDA  spot  size  diminish.  Far 
enough  downstream  the  shapes  of  the  profiles  will  reflect  the 
fluid  dynamics  and  not  be  influenced  by  the  LDA  dimensions,  as 
at  the  upstream  stations. 

Several  key  parameters  are  used  to  describe  the  jet  profile 
as  it  develops  downstream.  One  of  these  is  the  width  of  jet. 

(In  our  discussions  we  will  use  i>  1/2 <  the  full  width  at  half- 

maximum  velocity.)  There  are  many  references  to  the  development 

of  a  free  jet  into  the  same  fluid  (e.g.,  air-air).  Hinze20 

discusses  the  work  of  Forstall21  whose  measurements  extrapolate 

to  the  relation  $1/2  -  0.25  Z  for  a  free  jet  into  a  still  gas, 

and  the  work  of  Van  der  Hegge  Zijnen  who  obtained 

$1/2  =  0.16  (Z+a),22  where  a  is  the  virtual  origin  of  the  jet. 

Other  reported  measurements  are  $1/2  -  0.172  Z  23  and 

$1/2  =  0.186  Z  24  #  For  the  jet  studied  in  this  experiment, 

CO2  into  air,  the  results  shown  in  Figure  19  yield  the  best-fit 
relationship 

$  1/2  =  1.10  (Z/Dp  -  4.8)  (6) 

=  0.230  (Z  +  a) 


with 


a  =  -4.3  Dp. 

The  data  point  at  Z/Dp  =  104.7  was  not  used  since  it  was  deter¬ 
mined  that  at  this  location  the  central  jet  had  spread  so  much 
that  it  was  interacting  with  the  duct  boundary  layer. 
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AXIAL  LOCATION  (Z/Dp) 


Development 
Measured  at 
measurement 
the  jet  has 
the  duct. 


of  the  Full  Width  of  the  Central  Jet 
Half  the  Velocity  Maximum.  Mote  that  the 
at  Z./Dp  =  104.7  has  not  been  used  since 
spread  to  where  it  is  interacting  with 


Figure  19. 


Measurements  were  also  made  with  jet  flows  of  2  and 
16  kg/hr.  The  results  shown  in  Figure  19  encompass  all 
measurements.  The  profiles  at  2  kg/hr  (Figs.  20(a)  to  20(d)’  and 
16  kg/hr  [Figs.  21(a)  to  21(f)]  were  taken  at  a  limited  number  of 
stations.  They  were  used  for  comparison  with  the  profiles  of 
6  kg/hr  flow  to  investigate  systematic  variations  with  Reynolds 
number.  Figures  22(a)  to  (d)  are  overlays  of  the  averaqe  axial 
velocity  and  rms  axial  velocity  fluctuation  for  the  three  flow 
rates,  each  normalized  to  its  centerline  inlet  velocity,  at  the 
four  downstream  stations.  In  Figure  22(a),  the  inlet  profile 
shows  all  three  profiles  to  be  identical,  both  velocity  and  rms 
fluctuation.  In  Figures  22(b)  to  22(d)  the  profiles  of  the 
2  kg/hr  and  6  kg/hr  flows  remain  the  same;  the  16  kg/hr  flow 
profile  increasingly  deviates  from  the  other  two  the  farther 
downstream  the  measurement.  In  the  region  of  the  flow  field 
Z  >  4  cm,  all  the  profiles  have  the  same  shape;  the  flow  at 
16  kg/hr  does  not  attain  the  _same  proportion  of  the  maximum 
velocity.  The  shape  of  the  velocity  profile  is 


=  e-Z-29'z7DF-r-4-'2 


w 


(7) 


where  the  numerical  factor  in  the  exponent  was  selected  by 
matching  the  observed  width  data  of  Figure  19. 

Many  of  the  numerical  models  used  for  predicting  fluid  flow 
assume  isotropic  turbulence.  The  fully  developed  region 
illustrated  in  Figure  17  is  thought  to  be  the  region  charac¬ 
terized  by  isotropic  turbulence.  From  Figures  13(k)  and  21(f)  it 
can  be  seen  that  the  rms  fluctuations  have  not  become  isotropic 
at  Z/Dp  =  21  for  Re  =  31,000  and  Z/Dp  =  34  for  Re  =  33,000, 
considerably  farther  downstream  than  was  anticipated.  Comparing 
the  rms  fluctuations  of  the  V  component  in  Figures  13(h)  and 
1 3 ( i )  shows  that  the  radial  and  azimuthal  components  are  already 
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Figure  21  . 


(a)  -  (f)  Profiles  of  Velocity  Distribution 
Descriptors  3ased  Upon  the  First  Four  Moments. 
Measurements  are  for  the  centra  1- j at-c.oly  oond 
with  CO g  mass  flow  of  13  kg. 'hr. 
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igure  22.  Overlays  of  the  Normalized  Average  Axial-Velocity 
(Wr/W§0)  and  Normalized  RMS  Axial-Velocity 
Fluctuations  (w/w|0)  for  the  Flow  Conditions  of 
2,  5  and  15  kg/hr  CO2 •  W§0  is  the  average 

jet  axial  velocity  at  the  entrance  plane,  W 
the  average  axial  velocity  at  the  measurement  pla. 
and  w  is  the  average  axial  RMS  fluctuation  at  the 
measurement  plane.  The  profiles  were  measured  at 
Z/Dp:  (a)  o'.  4;  (b)  2.1;  (c)  3.4;  and  ii)  20.9. 
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identical  at  Z/Dp  =  10.5,  and  only  the  axial  component  is 
d  i f f eren t . 

There  are  several  possible  mechanisms  by  which  the  measure¬ 
ments  might  be  influenced.  First,  when  taking  LDA  measurements 
the  flow  must  be  seeded.  As  a  result,  the  seeded  flow  mixes  with 
unseeded  ambient  air.  In  flow  regions  where  the  flow  is  inter¬ 
mittent,  the  measurements  could  be  biased.  This  should  show  up 
in  a  time  analysis  of  the  velocity  record,  when  such  a  study  is 
made.  To  overcome  the  problem,  it  is  necessary  to  prepare  a 
stationary  seeded  ambient  atmosphere,  a  difficult  process  at 
best.  Second,  the  flow  is  being  measured  in  a  duct.  At  the 
farthest  downstream  stations  (Z/Dp  >  34)  the  jet  has  spread  to 
^ 1/2  >  0*30  Dq.  As  the  profile  broadens  it  will  interact  with 
the  duct  boundary  layer  flow  and  no  longer  be  a  simple  jet. 

Third,  the  mixing  and  spreading  of  the  jet  are  accompanied  by 
entrainment  of  the  surrounding  fluid.  This  fluid  comes  from  the 
duct  and  must  be  replenished  by  fluid  drawn  in  from  upstream  of 
the  bluff  body.  As  this  fluid  slowly  flows  past  the  centerbody, 
large  eddies  might  be  shed  which  could  influence  the  flow 
downstream.  Other  mechanisms  for  perturbing  the  flow  undoubtedly 
exist . 

The  measurements  shown  do  not  clearly  indicate  that 
perturbations  exist  within  the  measurement  region,  as  can  be  seen 
by  the  linear  development  of  the  profile  width  to  Z/Dp  =  105  as 

expected.  If  perturbation  does  occur,  the  on-axis  flow  should 

-F  -F 

show  the  effects  least.  Figure  23  is  a  pilot  of  W0/W00  vs  Z/Dc, 
the  centerline  velocity  normalized  to  the  issuing  velocity  as  a 
function  of  downstream  location,  for  the  three  flow  rates 
studied.  The  form  has  the  expected  shape,  a  flat  top  in  the 
potential  core  region,  a  transition  region,  and  then  a  region 
that  fails  off  as  Z“3.  In  the  region  Z/Dp  >  30,  the  centerline 
/elocity  decays  much  faster  than  anticipated. 
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The  16  kg/hr  results  show  that  the  apparent  Z-1  dependence 
occurs  upstream  of  the  2  and  6  kg/hr  velocity  decay.  This  is  in 
keeping  with  the  profiles  observed,  shown  in  Figures  22(a)  to  (d) 
and  discussed  earlier.  The  Z~1  dependence  region  is  usually 
described  by  the  relation^ 


Z +a 


This  is  a  first  order  approximation  to  the  relation 


and  the  simplification  is  only  valid  if  Z+a/Dp  >>  A.  The  apparent 
Z“1  dependence  observed  in  Figure  23  is  misleading.  Fitting  both 
forms  to  the  data  in  the  region  Z/Dp  >  12.6  yields  A  =  7.5  and 

a  =  -0.55  Dp  for  the  linear  approximation,  and  A  =  6.8  and 

a  =  -5.23  Dp  for  the  exponential  form  .  The  change  in  the  value 

of  A  is  small  but  the  shift  in  the  virtual  jet  origin,  a,  is 
large.  Both  curves  fit  the  data  within  the  expermiental  uncer¬ 
tainties  but  it  is  clear  that  the  criterion  Z+a/Dp  >>  A  is  not 
valid  until  Z/Dp  >  70,  beyond  the  region  of  measurement.  Thus, 
the  apparent  Z~1  behavior  is  misleading  and  the  linear  approxima¬ 
tion  should  not  be  used. 

Experimentally  determined  values  of  A  and  a  using  the 
velocity  decay  data  are  shown  in  Table  4.  Recall,  the  parameter 
a  was  also  determined  from  the  profile  width  data  by  Eq .  (6), 
a  =  -4.3  Dp. 
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TABLE  4 

PARAMETERS  DESCRIBING  CENTERLINE  VELOCITY  DECAY 


A 

a/DF 

Z 

L<  <  U 

Dp 

L  CJ 

REFERENCE 

5.9 

-0 . 5 

10 

50 

22 

5.9 

-3 

10 

50 

25 

5.4 

-7 

25 

90 

25 

7.5 

-6 . 5 

15 

50 

26 

6,8 

-5 . 2 

10 

80 

★ 

6.1 

-4.1 

10 

_ 

80 

*  ★ 

*  this  study,  Re  =  11,000*  and  Re  =  31,000 
**  this  study,  Re  =  83,000 


The  rms  velocity  fluctuations  are  plotted  in  Figure  24. 

They  exhibit  the  expected  behavior^7  rising  to  25  to  30% 
turbulence  intensity  at  Z/Dp  ~  30  and  holding  constant  farther 
downstream . 

The  velocity  profiles,  Figures  13,  20,  and  21,  were  combined 
with  the  CO2  concentration  measurements  of  Bradley, 19  (see 
Appendix  A)  made  under  the  same  flow  conditions.  A  numerical 
integration  scheme  was  used  to  obtain  the  total  mass  flow  within 
the  profile,  the  mass  flow  of  CO2  within  the  profile,  and  the 
momentum  flow.  The  results  are  shown  in  Figure  25,  where  the 
quantities  plotted  are  (momentum  f low )/( entrance  velocity) 2, 

(COg  mass  flow)/( entrant  mass  flow),  and  (total  mass 
f low )/( entrant  mass  flow).  The  particular  normal i zation  factor 
chosen  for  the  momentum  flow  has  no  special  significance.  (It  is 
equivalent  to  using  mass  flow.) 
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Figure  25.  Total  Mass  as  a  Fraction  of  Entrance  Mass  at  Each 
Downstream  Measurement  Station  (m  =  mc  +  mJ  .  The 
data  is  for  all  three  flow  rates  (2,  6,  anS  16  kg/hr 
C02).  The  normalized  CO2  mass  flow  and  normalized 
total  momentum  flow  at  each  station  is  shown  to 
demonstrate  the  consistency  of  the  results.  The 
dashed  lines  are  the  expected  values. 
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The  plot  of  the  total  mass  { normal i zed  to  the  entrant 
mass)  is  divided  into  two  regions.  In  the  region  from  the 
entrance  plane  to  Z/Dp  **  10  (plotted  with  the  symbol  x  in 
Figure  25)  the  entrained  mass  grows  linearly  with  the  downstream 
distance,  the  rate  of  growth  given  by 


m 


F 


m 


o 


(0.19  ±  0.03)  Z/Df  , 


and  m  =  m0  +  mp  is  the  total  mass  flow.  In  the  region  Z/Dp  >  20 
(plotted  with  the  symbol  o)  the  growth  rate  of  the  entrained  mass 
is  also  linear  but  the  rate  has  increased  to 


E 

o 


(0.25  ±  0.03)  Z/Dp  . 


Ricou  and  Spaldingl^  defined  a  constant,  C2 ,  describing  the 
growth  of  the  entrained  mass  which  they  found  to  be  independent 
of  the  density  of  the  jet 
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where  p0  is  the  jet  density  and  p s  is  the  density  of  the  gas  into 
which  the  jet  is  flowing.  Using  this  definition  the  experimental 
results  obtained  are 

C 2  =  0.24  ±  0.03  2/Dp  <  10 

=  0.32  ±  0.04  Z/Dp  >  20 

The  far-downstream  value  corresponds  to  the  result  obtained  by 
Ricou  and  SpaldinglS  and  Hill^-G.  But,  unlike  Hill  our  results 
show  a  linear  growth  region  near  the  entrance  plane.  Boguslawski 
and  Popiel,18  using  hot-wire  measurements  of  the  velocity 
profile,  also  observed  a  constant  value  of  C2  in  the  region 
0  <  Z/Dp  <  10,  but  their  measurement  yielded  C2  =  0.183, 
substantially  lower  than  obtained  here.  Their  study  was  for.  an 
air  jet  into  air.  Secondary  checks  and  momentum  flow  were  not 
presented  in  their  paper.  There  may  be  a  difference  when  a 
CO2  jet  is  used  in  this  initial  region  of  the  flow  development. 

In  the  far  downstream  region,  25  <  Z/Dp  <  500,  Ricou  and*  Spalding 
found  that  C2  was  independent  of  the  jet  gas,  but  this  may  not  be 
true  closer  to  the  entrance. 

It  should  be  noted  that  in  the  earlier  measurements  of  Ricou 
and  SpaldinglS  the  value  of  C2  depended  upon  the  Reynolds  number 
of  the  flow,  reaching  its  terminal  value  for  Re  >  25,000.  The 
measurements  being  reported  here  spanned  the  range  11,000  <  Re  < 
83,000.  No  consistent  dependence  upon  Re  was  observed. 

The  validity  of  the  measurements  was  checked  in  two  ways. 
First,  the  mass  flow  of  CO2  in  the  jet  should  be  constant  at  each 
downstream  station.  Figure  25  shows  that  at  every  downstream 
location  the  COg  mass  flow  in  the  jet  is  the  same  as  the  CO2  jet 
entrance  mass  flow  to  within  ±5%,  even  though  the  total  mass  has 
qrown  more  than  20  fold.  Second,  momentum  should  be  conserved 
for  the  jet.  The  momentum  divided  by  the  entrant  velocity 
squared  is  plotted  in  Figure  25.  The  normal ization  chosen  has  no 
special  significance  other  than  it  reduces  the  measurements  to  a 
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common  value,  for  the  nozzle  used,  independent  of  mass  flow. 

(For  more  generality  a  density  factor  should  also  be  used.)  The 
results  for  the  three  mass  flows  used,  are  shown  in  Figure  25. 
The  momentum  is  constant  to  within  t 10%  throughout  the 
measurement  region. 

The  consistency  checks  made  on  the  measurements  have  shown 
that  the  velocity  profiles  are  accrate  throughout  the 
measurement  region.  Therefore,  the  entrainment  measurements 
should  also  be  accurate.  They  indicate  that  the  entrainment 
changes  notedly  when  the  jet  is  in  the  transition  region  of  its 
development  (see  Fig.  17). 
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SECTION  4 
ANNULAR  FLOW 


The  annular  jet  flowfield  provides  the  numerical  models  with 
a  substantive  test  for  evaluating  their  effectiveness.  A  recir¬ 
culation  vortex  is  established  in  the  wake  of  the  bluff  body 
(Figure  2)  and  the  reverse  flow  on  centerline  opposes  any  flow 
from  the  central  jet.  The  flowfield  was  measured  with  an  annular 
flow  rate  of  2  kg/s  of  air  and  with  central  jet  flows  of  0  and 
6  kg/hr  of  COg.  The  initial  results  reported  are  for  the  case 
with  no  central  jet  flow.  The  addition  of  the  central  jet  pro¬ 
duced  a  localized  change  in  the  flowfield  which  did  not  extend 
far  downstream.  A  detailed  study  of  the  perturbation  is  reported 
next . 

4.1  EXPERIMENTAL  PROCEDURE 

The  flowfield  was  measured  using  the  LDA  instrument  pre¬ 
viously  described.  The  annular  jet  and  central  jet  are  seeded 
independently,  and  care  must  be  taken  to  match  the  seeding 
densities  as  closely  as  possible. 

The  experimental  procedure  for  measuring  the  flowfield 
consists  of:  first,  scanning  the  flowfield  in  the  X  direction 
at  Y=0  (see  Figure  1);  second,  taking  measurements  in  the  Y 
direction  at  X=0  to  the  limit  of  the  window  (at  Y=20  mm);  and, 
last,  scanning  the  flowfield  in  the  X  direction  at  Y=20  mm 
(Fig.  26).  The  first  scan  provides  the  axial  and  azimuthal 
velocity  components.  The  scans  are  taken  from  one  side  of 
centerline,  through  centerline  to  the  window  on  the  other  side. 
The  scan  region  around  centerline  permits  checking  of  the  flow 
symmetry.  The  second  scan  gives  the  axial  and  radial  velocity 
components.  The  third  scan  determines  the  axial  velocity  and  a 
combination  of  the  radial  and  azimuthal  components  (Fig.  26). 
Assuming  that  the  average  flowfield  is  stationary  and  that  the 
average  azimuthal  velocities  are  known  from  the  first  scan,  the 
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figure  26. 


On-Axis  View  of  APL  Combustor  Showing 
Profiles  Taken  at  the  Axial  Stations, 
measures  v  and  w  components • 


the  Three  Sea 
The 
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radial  component  can  be  obtained  from  the  combined  measurement  by 
the  decomposition  given  in  Eqs.  (3a)  to  (8c). 


R  = 

/x2  +  4 

(3a) 

a  = 

tan-1 ( 2/X) 

(3b) 

VR  = 

(sin  9)“^  (V  -  V  cos  9). 

9 

(3c) 

where  (R,  9)  are  the  cylindrical  coordinates  of  the  measurement 
location  at  (X,2)  in  cartesian  coordinates, 

Vr,  Vg  are  the  radial  and  azimuthal  velocity  components, 
and  V  is  the  velocity  component  in  the  Y  direction. 

4.2  RESULTS 

The  measurements  taken  at  Y=0  are  shown  in  Figures  27(a) 
to  (n) .  The  measurements  were  made  at  locations  0  <  Z/Dc  <  4.29, 
where  Dc  is  the  centerbody  diameter,  140  mm.  The  edge  of  the 

9 

centerbody  occurred  at  r/Rp  =  29.3.  The  relatively  flat  entrance 
profile  of  the  annular  jet  is  seen  in  Figure  27(a).  The  rms 
velocity  fluctuations  are  uniform  until  the  centerbody  boundary 
layer  is  entered,  and  then  they  increase  to  four  times  the  value 
in  the  annulus.  As  the  measurements  proceed  downstream,  the 
recirculation  zone  increases  in  strength,  the  rms  fluctuations  in 
the  recirculation  zone  increase,  and  the  fluctuations  in  the  shear 
layer  diminish.  This  process  continues  until  at  Z/Dc  =  0.29  the 
shear  layer  is  not  distinguishable  from  the  recirculating  flow  on 
the  basis  of  the  fluctuating  component. 

In  the  region  Z/Dc  *0.86  the  fluctuating  components  have 
reached  their  maximum.  Further  downstream  the  axial  fluctuation 
falls  faster  than  the  azimuthal  fluctuation  and  this  remains  the 
case  to  the  end  of  the  first  window's  observations  at 
Z/Dc  =  1.43.  In  the  second  window  the  velocity  profile  is 
approaching  a  flat-top  profile  and  the  fluctuating  components  are 
equal  for  the  two  directions  measured. 
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The  other  two  scans  were  combined  with  Figure  27  to  provide 
the  radial  velocity  profiles  shown  in  Figure  23,  according  to  the 
mathematical  manipulation  of  Sqs.  (3a)  to  (3c).  These  profiles 
should  be  used  cautiously.  Once  the  measurement  location  is 
beyond  r/Rp  ~  10  the  (sin  6)-^  factor  quickly  becomes  large  and 
can  combine  with  small  errors  in  V  to  severely  distort  the 
profile.  The  velocities  reported  are  small  and  the  uncertainties 
are  comparatively  large.  The  profiles  should  oe  used  to  indicate 
trends.  They  do  display  the  expected  inversion  at 
*0.29  <  Z/Dc  <  0.43  indicating  planes  on  either  side  of  the  vortex 
center . 

The  profiles  shown  in  Figures  27(a)  to  (n)  were  input  to  the 
mass  and  momentum  computation  program  discussed  in  the  section  on 
free  jets.  The  results  are  shown  in  Figure  29.  Also  plotted  are 
the  results  using  the  coaxial  flow  profiles,  the  addition  of  a 
6  kg/hr  central  jet  having  no  observable  impact  on  the  2  kg/s 
annular  jet  contribution.  The  computed  mass  flow  is  2  kg/s  ;3% 
over  the  entire  range  of  axial  stations  0  <  Z/Dc  <  4.29.  The 
axial  momentum  flow  starts  at  the  value  expected  for  a  flat-top 
velocity  profile  in  the  annulus.  As  the  flow  in  the 
recirculation  zone  builds  up,  the  momentum  flow  diminishes.  The 
decrease  continues  past  the  air  stagnation  position,  until  the 
flow  in  the  region  blocked  by  the  centerbody  has  .built  up  to 
match  the  decreased  velocity  around  the  periphery.  At  this  point 
the  flow  appears  to  be  almost  a  flat-top  velocity  profile  across 
the  duct,  and  the  momentum  bears  this  out.  As  the  flow  proceeds 
downstream,  it  is  expected  to  become  fully  developed  pipe  flow 
across  the  duct.  The  anticipated  momentum  flow  is  shown  on 
Figure  29.  Typically  a  length  of  Z/Dc  of  50~100  is  required  to 
achieve  this  condition  so  it  is  not  observed  in  the  modest 
downstream  region  observable  (Z/Dc)max  =  4.29  although  some 
recovery  is  observed. 

The  mass  and  momentum  program  provides  an  output  of  the 
integrated  mass  between  centerline  and  each  radial  step  out  to 
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the  duct  wall.  These  data  are  then  interpolated  to  determine  the 
radial  locations  which  encompassed  specific  mass  flow  fractions. 
3y  connecting  these  locations  of  constant  mass  flow  fractions  the 
streamlines  for  the  flow  are  determined.  These  are  shown  in 
Figure  30.  The  region  of  the  centerbody  0  <  r/Dc  <  0.5  is  shown 
shaded  on  the  axis;  0.5  <  r/Dc  <  0.91  is  the  annulus.  The 
recirculation  region  contours  were  evolved  with  substantial 
artistic  license  in  the  regions  of  severe  curvature  where  data 
points  are  sparse.  It  should  be  remembered  that  the  data  points 
are  measured  along  planar  cuts  through  the  flow.  Consequently, 
it  is  not  possible  to  follow  a  streamline  to  trace  it  out.  Only 
one  planar  cut  showed  a  reverse  streamline  of  -0.065  and  this  is 
shown  as  a  point  in  Figure  30.  But,  from  the  radial  velocity 
data  shown  in  Figure  28,  the  center  of  the  vortex  should  be 
upstream  of  this  point.  Unfortunately,  the  planes  of  measurement 
are  too  far  apart  to  locate  the  -0.065  contour  and  the  vortex 
center  cannot  be  determined  any  more  precisely  than 
0.36  <  Z/Dc  <  0.43,  0.34  <  r/Dc  <  0.38.  • 

Radial  scans  along  the  X  axis  (Y=0)  were  made  with  the  flow 
conditions;  central  jet  6  kg/hr  CC>2?  and  annular  jet  2  kg/s 
air.  These  were  taken  at  axial  positions  0.01  <  Z/Dc  <  0.57,  at 
which  point  it  was  determined  that  there  was  no  observable 
difference  between  the  measurements  with  and  without  the  central 
jet.  The  influence  of  the  central  jet  is  restricted  to  a  region 
near  the  centerline  and  this  was  examined  in  more  detail,  at  each 
axial  station,  by  scanning  in  the  Y  direction  at  X=0 .  These 
measurements  are  also  shown  in  Figure  31.  At  the  entrance  plane 
and  for  a  short  distance  downstream  0  <  Z/Dc  <  0.04  the  central 
jet  has  similar  characteristics  to  those  measured  for  the  free 
jet.  The  jet  profile  decays  much  more  rapidly  due  to  the  reverse 
flow  from  the  annulus,  and  the  velocity  fluctuations  increase 
more  dramatically.  At  the  measurement  plane  Z./Dc  =  0.29  the  jet 
is  barely  discernable,  but  the  axi3l  velocity  fluctuations 
dramatically  give  evidence  to  the  presence  of  the  central  jet. 
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Compare  Figures  31(e.l)  and  27(e).  At  Z/Dc  =  0.43  and  beyond, 
the  Clow  properties  have  returned  to  those  measured  with  no 
central  jet. 

Centerline  profiles  were  made  for  the  average  velocity  under 
the  Clow  conditions:  annular  jet  Clow  of  2  kg/s  air  and  central 
jet  flow  of  0,  2,  6,  16  kg/hr  CO2.  They  are  shown  in  Figure  32. 
In  all  cases  measured,  reverse  flow  is  observed  in  the  average 
velocity.  The  maximum  reverse  velocity,  about  41%  of  the  annular 
inlet  plane  velocity,  is  the  same  as  was  measured  previously . 26 
The  rms  fluctuations  are  shown  in  Figure  33. 

The  locations  of  the  central  jet  stagnation  position  were 
plotted  against  the  ratio  of  the  central  jet  inlet  velocity  to 
annular  jet  inlet  velocity.  The  result  is  shown  in  Figure  34. 
Also  plotted  are  the  results  of  the  previous  1-D  LDA  measurement, 
using  a  different  nozzle. 26  it  can  be  seen  that  with  the  nor¬ 
malization  used  the  curves  are  almost  identical.  Thus  it  appears 
that  the  location  of  the  fuel  stagnation  point  scales  with  the 
ratio  and  not  with  the  mass  or  momentum  of  the  central  jet. 

and  wj  are  the  central  jet  and  annular  jet  entrance  plane 
velocities. 
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SECTION  5 
CONCLUSIONS 


A  two-dimensional  Laser  Doppler  Anemometer  has  been  designed 
and  constructed  for  use  in  a  ducted  coaxial  flow  system  having  a 
centerbody.  The  LDA  has  been  used  to  take  measurements  in  the 
flowfield  in  three  flow  conditions:  central- jet-only ;  annular- 
jet-only;  and  coaxial  flow.  The  measurements  have  provided  the 
information  needed  to  compute  the  mean  value  of  the  velocity  com¬ 
ponent  and  the  fluctuating  component.  In  addition  higher  moments 
of  the  distributions  have  been  computed. 

The  parameters  describing  the  central- j et-only  measurements 
have  been  compared  with  the  standard  descriptions  obtained  by 
past  investigators  who  have  used  photography,  hot-wire,  and  other 
study  techniques.  The  central- jet-only  measurements  have  also 
been  used  with  concentration  measurements  made  by  intrusive 
probes  under  the  same  operating  conditions,  to  compute  the 
entrainment  of  the  surrounding  fluid  into  the  jet. 

The  jet  half-maximum-velocity  width  has  the  usual  develop¬ 
ment  shape,  but  the  asymptotic  linear  growth  for  this  CO2  into 
air  jet  is  at  a  much  higher  rate  than  observed  for  the  air  jet 
into  air  studies,  $i/2/(z+a)  is  0.262  vs  0.17.  The  entrainment 
factor  C2  is  constant  in  these  studies  over  the  range 
0  <  Z/Dp  <  10,  not  showing  the  growth  region  observed  by  Hill, 
and  the  value  observed  0.24  is  larger  than  that  obtained  by 
Boguslawski  and  Popiel  (0.18)  for  an  air  jet  into  air.  Also,  no 
marked  dependence  of  C2  on  Reynolds  number  was  observed  over  the 
range  studied  11,000  <  Re  <  83,000,  in  contrast  to  the  earlier 
results  of  Ricou  and  Spalding,  although  their  measurements  were 
in  a  region  much  farther  downstream.  Beyond  Z/Dp  ~  10  the  value 
of  C2  increases  to  0.32  matching  the  value  obtained  by  Ricou  and 
Spalding.  Consistency  checks  made  on  the  central  jet  have  shown 
that  the  mass  flow  of  CO2  is  constant  downstream,  even  in  the 
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presence  of  a  20  fold  increase  in  mass  by  the  air  entrained  and, 
also,  the  jet  momentum  is  conserved.  Both  these  results  give 
confidence  in  the  measured  values. 

The  annular- jet-only  flowfield  has  been  mapped  extensively 
in  the  region  of  the  wake  of  the  centerbody  and  downstream  from 
it.  The  recirculation  vortex  has  been  studied  and  the  center 
located  in  the  region  0.34  <  r/Dc  <  0.38  and  0.36  <  Z/Dc  <  0.43, 
and  the  maximum  strength  streamline  observed  had  the  value 
-0.065.  The  recirculation  zone  extends  downstream,  on  centerline 
for  a  distance  Z/Dc  ~  0.97.  To  test  the  consistency  of  the 
measurements  the  mass  flow  was  computed  for  each  velocity  profile 
measured  and  at  every  station  the  value  obtained  was  the  inlet 
mass  ±3%.  The  axial  momentum  flow  of  the  annular  jet  started  at 
the  value  of  the  annular  flow  flat-top  entrance  profile.  It 
decreased  downstream  until,  at  a  position  well  past  the 
stagnation  point,  it  approached  the  value  describing  a  flat-top 
velocity  profile  across  the  duct.  After  this  point  it  began  to 
increase,  but  in  the  relatively  short  length  of  duct  observable, 
Z/Dc  =  4.29  maximum,  it  did  not  have  time  to  reach  the  value 
associated  with  fully  developed  pipe  flow. 

The  coaxial  flow  profile  measurements  traced  the  decay  of 
the  central  jet  due  to  the  action  of  the  reverse  jet  set  up  by 
the  annular  flow  recirculation  vortex.  As  the  velocity  change 
resulting  from  the  central  jet  diminishes,  the  rms  fluctuation 
increases  strongly,  indicating  the  great  amount  of  mixing 
occurring  in  the  region  near  the  central  jet  stagnation  location. 
In  combustion  studies,  a  flame  ball  is  usually  observed  in  this 
region  giving  visual  evidence  to  the  mixing.  Centerline  profiles 
were  measured  for  three  central  jet  flow  rates,  and  the  results 
obtained  showed  good  agreement  with  previous  measurements 
obtained  with  a  one-dimensional  LDA. 
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APPENDIX  A 


C02  CONCENTRATION  MEASUREMENT  PROFILES 
FOR  CENTRAL  JET  FLOWS  INTO  AIR 


APPENDIX  A 


The  CO2  concentration,  measurements  of  Bradley,^  in  the 
region  Z/De  >  20,  have  been  fit  to  functional  dependence  of  the 
same  form  as  the  LDA  measurements.  The  result  is 


where  Ccq?  is  the  fractional  concentration  of  CO2  in  the  gas 
sampled  at  location  (r,  Z)  [both  in  mm]  . 

In  the  same  region  the  LDA  was  shown  to  have  the  rela¬ 
tionship  (for  2  and  6  kg/hr  of  CO2) 

"6-8  Z/D  5.2  ~2*29  (z/d/-  4.3) 2 

W  ( r ,  Z)  =  ( 1  -  e  )  e  '  j.(A.2) 

The  general  fit  to  the  data  is  quite  good.  There  are  some  points 
though  that  should  be  noted.  First,  the  numerical  coefficients 
in  the  two  equations  are  not  the  same  as  would  be  expected  for 
fully  developed  turbulent  flow.  Second,  the  'ocation  of  the  vir¬ 
tual  origin  is  different  in  each  exponential  expression  contrary 
to  what  was  expected. 20  The  discrepancies  do  -oh  appear  to  be 
unreconcilable .  They  may  be  due  to  experimental  arror  and/or  to 
the  fact  that  the  data  used  was  not  from  a  reqion  far  enouqh 
downstream  to  represent  the  truly  self-similar  region  of  the 
fluid  flow.  Further  study  will  be  required  to  resolve  the  issue. 
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